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Abstract 

Nitrogen doped graphene and carbon nanotubes are popularly in focus as metal-free electro-catalysts for 
oxygen reduction reactions (ORR) central to fuel-cells. N doped CNTs have been also reported to chemisorb 
mutually, promising a route to their robust pre-determined assembly into devices and mechanical reinforcements. 

We propose from first-principles a common mechanistic understanding of these two aspects pointing further 
to a generic chemical activation of carbon atoms due to substitution by nitrogen in experimentally observed 
conhgurations. Wannier-function based orbital resolved study of mechanisms suggests increase in C-N bond- 
orders in attempt to retain 7r-conjugation among carbon atoms, causing mechanical stress and loss of charge 
neutrality of nitrogen and carbon atoms, which remedially facilitate chemical activation of N coordinated C 
atoms, enhancing sharply with increasing coordination to N and proximity to zigzag edges. Activated C atoms 
facilitate covalent adsorption of radicals in general, diradicals like O 2 relevant to ORR, and also other similarly 
activated C atoms leading to self-assembly of graphenic nano-structures, while remaining inert to ordinary 
graphenic C atoms. 

Disruption of 7r-conjugation due to substitution by p— or n— type dopants in graphene causes localization of 
2pz electrons, which are a rich source of exotic physical as well as chemical functionalities unavialable in undoped 
graphene. Accordingly, boron (B) and nitrogen (N) doped graphene (GF) and carbon nanotubes(CNT) have 
been under rigorous scrutiny during the last decade or so, resulting into a great multitude of proposals for novel 
applications ranging from nano-electronics to catalysis. Electro-catalytic functionality of doped GFs and GNTs 
have been under investigation (Hiiaisi) for more than a decade, primarily in pursuit of an efficient metal-free 
catalyst to replace the expensive platinum based ones used in fuel-cells so far. In particular, N doped GNTs 
and GFsdi 0 n H 0 m uni El E|), and more recently their composites (pRl [T4l [T5l [T6l ITT)) . have been explored 
exhaustively in the last few years as catalysts for oxygen reduction reactions (ORR) promisimg removal of platinum 
from the cathode in fuel-cells. On the other hand, N doping has also been suggested (jlS]) to induce covalent mutual 
adsorption among GNTs through C-G cross-links, as has indeed been reported(|T9]), promising a new route to 
their controlled assembly into devices and stronger composites, which in itself is an important open problem. 
In the present work we show mechanistically that indeed these two aspects are just different facets of a more 
inclusive general scenario of chemical activation of G atoms arising with experimentally observed configurations of 
substitutional doping by N in sp‘^ hybridized carbon (C) based graphenic nano-structures (GN). 

First principles studies ([^ 1^ [2^ 1^ of ORR catalysis so far, primarily of isolated mono-substitution by N, 
suggest adsorption of atomic, molecular as well as active anionic oxygen (O) on C atoms in the vicinity of quaternary 
N atoms, presumably driven by the effective positive charge of G atoms due to electron extraction by their N 
neighbours, which however does not explain the observed mutual adsorption of N doped GNTs. Interestingly, recent 
direct observations suggest substitution by N in GF to occur predominantly at next-nearest neighbouring 

(N-Nn) sites, which hinders one of the three degenerate 7r-conjugation configurations more than the other two, unlike 
in case of sparse mono-substitution where all three of them are equally hindered. This implored us to ask whether 
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Figure 1: (a) Energetics of substitutional doping, (b-e) Variation in bond lengths compared to corresponding 
equilibrium values. (f,g) Contour plot of the degree of non-co-planarity measured as the vertical distance of an 
atom from the plane passing through its nearest neighbours, (h) Histogram of the degree of non-co-planarity for 
different substitution configurations. 

multiple substitution at N-Nn sites can open up a new class of mechanisms supporting ORR catalysis as well as 
GN-GN covalent adsorption on the same footing. Additionally, as the N-Ns sites belong to one of the two sub¬ 
lattices of graphene, which support inter-sub-lattice spin-separation upon physical or chemical disparity between 
the two sub-lattices, the resultant magnetism needs to be linked to other functionalities arising out of substitution 
at N-Nn sites. Notably, inter-sublattice spin-separation, which manifests itself as the nearest neighbour anti-ferro¬ 
magnetic order, and can be rationalized as an effect of on-site Goulomb repulsion between electrons with opposite 
spins, is naturally expected to contest 7r-conjugation, although modestly. Disparity between sub-lattices due to 
zigzag edges, which is well known to consolidate nearest neighbour anti-ferro-magnetism in GF, is thus expected 
to considerably influence properties of GNs with substitution by N at N-Nn sites. 

Through unambiguous quantitative estimates of spin resolved bond orders and atomic charges obtained based 
on spatially localized Wannier functions(WF), we find the onset of mechanical strain and loss of charge-neutrality 
of N and G atoms due to evolution of the G-N bond-orders in attempt to maximally retain 7r-conjugation among 
G atoms upon substitution by N, to not only govern the energetics of substitution but also to favour passivation 
of N coordinated C atoms as an overall remedy. Such activated C atoms facilitate covalent adsorption of radicals, 
di-radicals like O 2 , as well as similarly activated G atoms leading to cross-linking of N doped GNs, although remain 
inert to ordinary graphenic G atoms, which is advantageous since it would allow spatially resolved self-assembly of 
N doped GNs. Adsorption of O 2 as diradical leaves the mono-coordinated O atom charged and reactive as required 
for ORR. Activation is found to increase sharply with increasing coordination to N and proximity to zigzag edges. 
Absolute magnetization is observed to alter significantly in events of adsorption, since it restores 7r-conjugation to a 
variable degree, which in turn modifies nearest neighbour anti-ferromagnetism, given the anti-correlation between 
the two. In addition to proposing a new paradigm of activation of graphenic G, this work also highlights the 
effectiveness of WFs in shaping precise orbital resolved understanding of mechanisms of chemical processes. 

Equilibrium configurations and relevant energetics of substitutional doping are obtained within the non-empirical 
framework of density functional theory (DFT) based on local mean-field approximation of the many-electron 
exchange-correlation contributions to total energy. We use the Quantum Espresso ([27]) code, which expands wave- 
functions in the basis of plane waves and allow ultras 0 ft pseudopotentials for valence electrons. For exchange- 
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correlation we use a gradient corrected Perdew-Burke-Ernzerhof (PBE)(j28j) functional. We considered a 5x5 GE 
super-cell to understand the energetics of substitution by N and the proposed activation of C atoms thereby. To 
study the effect of zigzag edges and use an available scheme for unbiased partitioning of total charge density in 
isolated systems, we also consider an isolated GE segment (EigQ a)) made of 84 atoms. Total energies and con¬ 
figurations are converged with plane-wave cutoffs over 800 eV, k-mesh up to 7x7x1 for the super-cell, and forces 
less than 10“^ Rydberg/Bohr using the Broyden-Eletcher-Goldfarb-Shanno (BEGS) ([30]) scheme for total energy 
minimization. To account for the attractive dispersion interactions, semi-empirical Grimme approximation ([3T]) is 
used only in cases where chemisorption is not conclusive with PBE. 

Unambiguous estimates of bond-orders, atomic charges and sub-shell filling, are obtained through construction 
of localized WEs based on joint-diagonalization of non-commuting first moment matrices (EMM). In one dimension, 
WEs with maximum localization([34|) are naturally the eigen-functions of the EMM which is essentially the position 
operator expanded in the occupied subspace. Eor periodic systems the EMM is obtained in terms of geometric 
phases of Bloch states evolving over the full Brillouin zone. Eor isolated systems the EMMs can be directly 
calculated from the occupied Kohn-Sham (KS) eigen-states as 

( 1 ) 


for occupied KS states in the valence band. Diagonalizing a EMM (say X as evaluated above) yields WEs with 
maximum localization in the corresponding {x) direction. However, owing to the inherent non-commuting nature 
of the position operators expanded within the finite occupied subspace, it is not possible to obtain WEs with 
maximum localization simultaneously in all three directions. Nevertheless, approximate joint diagonalization of 
the three EMMs renders a set of highly localized WEs, which, even though strictly not unique, constitute an 
unambiguous orbital representation of the occupied sub-space. Their non-uniqueness arise only from the sequence 
in which the three EMMs are considered in course of the iterative approximate joint diagonalization process, and 
in effect has negligible bearing on the WEs rendered. Wannier centres (WGs), which are centre of masses of WEs, 
are readily available as the approximate eigen-values of the three EMMs, without having to explicitly construct the 
WEs. WGs, each representing one electron for each spin, provide a unique dot structure map for valence electrons 
through out the system. In case of partially occupied bands, the map is obtained by weighted sum of WGs for 
different number of KS states: 
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where /jv is the occupancy of the A^-th KS state, such that 
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where Ne is the total number of valence electrons in the isolated system. 

WGs can be identified in two categories - (I) associated with atoms, and (2) representing bonds. Single and 
double bonds are represented respectively by one and two WGs between two atoms. Eor visualization purpose, WGs 
in close proximity, like the ones representing double bonds, have been fused into one WG representing two electrons. 
WGs are represented by spheres, whose radii are scaled linearly by the number of electrons they represent, which 
is thus between 0 to 2. The depth of shading(in gray-scale) of the planar projection of these spheres are also 
linearly scaled in the same range. Number of electrons associated with an atom is implied by the WGs exclusively 
associated with it, representing unpaired and lone-paired electron(s) if any, and half the order of all the bonds 
made by the atom according to the WGs representing covalent bonds. WEs and WGs are constructed using an 
in-house implementation which takes the KS eigen states as input. The reason to use this method is that it does not 
depend on any reference template of orbitals, which is required for construction of localised WEs ([35l[36|) in periodic 
systems, where it also biases the resultant WEs towards the set of orbitals specified in the template. However, 
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Figure 2: (a-k) Planar projection of WCs implying order of bonds, charge state, spin polarization and subshell¬ 
filling of atoms. Adjacent WCs have been merged into one WC located at their center of mass. Larger and darker 
circles each denote 2 electrons per spin, while the smaller gray circles represent 1 electron per spin. (1) Spin density 
{pspini — Pspin2) with substitution at three mutually N-Nn sites enclosing a carbon atom. 


as the method can not be directly extended to periodic systems, estimation of bond-orders and atomic charges 
are made based on WCs calculated in hydrogen passivated isolated GF segments, and the implied mechanisms 
are argued to be plausible in their corresponding periodic super-cells based on the similarities in bond-lengths, 
bond-angles and the degrees of non-co-planarity of atoms observed in the periodic super-cells and their isolated 
GF counterparts. 

Energy required for substitution by N in the periodic super-cell, calculated as {Edoped — Eundopped) + mEc — 
ttiEn), where E is total energy and m is the number of C atoms substituted, is plotted in Figj^a). Gonsistent 
with experimental observations m m, Fig§ a) implies substitution at nearest neighbouring (Nn) sites to be 
particularly unfavourable. The deviation in equilibrium bond-lengths plotted in Fig{^b-e) from that of the res¬ 
onating graphenic G-G bonds and single (order 1) C-N and N-N bonds obtained with same parameters, suggests 
strain in the G-N bonds to be possibly responsible for difference in energetics of substitution. To translate these 
deviations into strain, bond orders have to be estimated, since reduction in length of a bond can also happen with 
increase in its order.Fig J^e) suggests complete retention of 7r-conjugation among the C atoms with substitution at 
Nn sites supporting C-N and N-N bonds or order 1, as readily understood by considering substitution at the two 
sites enclosed by the dashed rectangle in Figj^a). The C atoms would thus strongly favour co-planarity among 
themselves, while the N atoms, owing to their lone pairs, would strongly favour non-co-planarity with their nearest 
neighbours. We indeed observed this in the periodic super-cell, as evident in Figj^f), which maps the degree of 
non-coplanarity of an atom as its vertical distance from the plane passing through its three nearest neighbours. 
Thus with substitution at Nn sites, the C-N bonds in the periodic supercell are also likely to be of order 1, implying 
that those C-N bonds are indeed longitudinally strained as their observed reduction in length is not due to any 
increase in their order. The stark contrast between C and N atoms in their preference for non-coplanarity with 
their nearest neighbours also suggests non-nominal dihedral strain. 

With substitution at non-Nn sites, completely unaltered retention of 7r-conjugation is impossible, resulting into 
ease in accommodation of C-N single bonds, which are typically about 4% longer than the graphenic C-C bonds. 


4 




























Figure 3: (a) Energetics of adsorption of CH 3 on different C atoms, (b) Interaction potential of CH 3 as a function 
of vertical distance of its C atom from the C atoms in periodic GF supercells specified in the figure, (c) Equilibrium 
configurations of adsorbed CH 3 on Co, C’ and C* in periodic GF supercell. 

This is evident in Figj^d) which shows only two of the six C-N bonds to have significant reduction in length upon 
bi-substitution at N-Nn sites. Furthermore, Fig§c,h, i) suggests increase in order of C-N bonds with substitution 
at non-Nn sites, despite the consequent loss of charge neutrality of N atoms due to depletion from their lone pairs. 
This in effect imply that retention of 7 r-conjugation among C atoms is the dominant mechanism which can force 
locally unfavourable transfer of charges to make way for global lowering of total energy. Depletion from lone pairs 
of N would also imply reduced preference for non-coplanarity of N atoms with their nearest neighbours, which is 
indeed evident from Figj^h) as it suggests lesser degree of non-coplanarity of N atoms in periodic supercells with 
non-Nn substitution. The corresponding C-N bonds in the periodic supercell are thus expected to be of order i 1 
which would account for their reduced length (Fig{^b,d,e)), implying that they are effectively unstrained and thus 
energetically favourable. 

Hindered 7 r-conjugation is expected to allow the nearest neighbour anti-ferromagnetic correlation to consolidate, 
which is further expected to reinforce due to electron localization induced by the zigzag edges in the vicinity. In the 
isolated segment spin-separation is observed (FigJ^b-c,f-i)) with all possible non-Nn substitution, whereas in the 
periodic supercell spin separation consolidates (FigJ^l)) only with tri-substitution at mutually N-Nn sites enclosing 
a common nearest neighbouring C atom(C* in FigJ^. This again reiterates that the C atoms would always like 
to maximize 7r-conjugation among themselves, which can be resisted effectively by the nearest neighbour anti- 
ferromagntic order only when all the three 7r-conjugation configurations are equally compromised, as happens with 
enclosing tri-substitution, or in the presence of magnetic impurity or zigzag edges. Within the region of spin- 
separation each C atom has in effect an unpaired electron centred on them (FigJ^b-c,f-i)), indicating their partial 
sp^ hybridization, which is expected to help cooperatively in allowing the local non-coplanarity preferred by the 
N atoms. However, as the N atoms themselves have reduced preference for non-coplanarity due to depletion from 
their lone pairs, an effectively planar equilibrium configuration is attained with mild systematic corrugation due to 
the partial sp^ hybridization of anti-ferromagnetically correlated C atoms, as evident in FigQg) for the periodic 
supercell with enclosing tri-substitution. 

Energetics of adsorption of Me (CH 3 ) radical on C atoms shown in FigJ^a) suggests adsorption on the N 
coordinated C atom (Co) due to enclosing tri-substitution to be prominently most favourable. As evident from 
Fig§b, c), mono-substitution by N in the isolated segment leads to spin-separation and a C-N double bond for 
the majority spin. Interestingly, passivation of a C (C* in FigJ^ atom nearest to the N atom, exemplified here 
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Figure 4: Adsorption of NO 2 (a), NO (b) and O 2 (c) on Cq. (d-e) WCs corresponding to (a),(b) and (c) respectively, 
depicted by smaller (larger) gray (black) circles denoting 1(2) electrons per spin, (g-h) WFs representing the Cq-Co 
cr-bond (g) and an inter-segement tt-tt bond (h). (i) Inter-CNT Cq-Co and adjoining C-C bonds. N atoms are 
marked by black dots in (g-i). 


through adsorption of Me, not only restores complete 7r-conjugation by quenching spin-separation, but also reduces 
the order of the C-N double bond (Figj^c)) to 1, which restores charge neutrality of the N atom and relieves 
the tensile stresses in the adjacent grapheme C-C bonds. In case of substitution at Nn sites, which retains tt 
conjugation among C atoms intact, passivation of any C atom beyond the nearest neighbourhood of N atoms is 
same as passivation of a C atom in an undoped segment as evident in FigJ^a). Passivation of C atoms (C” in 
FigJ^ next to an N atom is more favourable on account of cooperative non-planarity of the N and the passivated C 
atom, although it hinders tt conjugation locally and consolidates nearest-neighbour antiferromagnetic order in the 
vicinity. With the experimentally observed substitution at N-Nn sites, the C atom (C’ in FigJ^ nearest to both 
the N atoms, neither retains charge neutrality nor completes sub-shell filling. Interestingly, both the shortcomings 
are restored upon passivation of C’, which also makes the C atom next to it charge neutral, implying a higher level 
of activation of C’ than an ordinary grapheme C atom, as also evident in FigJ^a). 

With enclosing tri-substitution, order of Cq-N bonds (Figj^h,i)) increases to around 1.5 on the average, implied 
by the three C-N double bonds out of six total C-N bonds accounting for the two spins separately. This is consistent 
with the shortest C-N bond lengths observed in the periodic super-cell among the different possible N substitution 
configurations, implying high degree of mechanical stress in the adjacent C-C bonds. While strong depletion 
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Figure 5: (a) Distance between the lower O of O 2 , and the active C atoms as shown in the upper panel, in course of 
structural relaxation using the BFGS (|3Q|) scheme. DFT + Dispersion indicates use of Grimme’s approximation (|3T]) 
along with PBE exchange-correlation to account for the dispersion interation. (b-c) Free energy diagram for 
complete reduction of the diatomic oxygen adsorbed on the active G atoms shown in the inset. These free diagrams 
are to be compared with free diagram (|4Q|) for ORR on a G atom next to a quaternary N atom, which shows an 
overpotential of about 0.7 V. 

from lone-pairs of N atoms leaves them positively charged and helps the nearest neighbour anti-ferromagnetism 
to consolidate. Go would have excess charge centered on it owing to the increased order of the Gq-N bonds. This 
excess charge on Co, being in excess to what it requires for complete sub-shell feeling, would increase total energy 
and should have low ionization potential. Conveniently, as evident in (Figj^j,k)), complete sp^ hybridization 
of Co, as possible upon passivation, would force all Co-N single bonds to be of order 1, which will restore charge 
neutrality of of N atoms, besides reducing tensile stresses in the surrounding C-C bonds and facilitating cooperative 
non-planarity of Co and N atoms. Co would thus be prone to passivation, and thereby, non-trivially active in 
terms of adsorption of radicals, as evident from the adsorption energies plotted in Figj^a). Interaction potential 
(^(Me+GF) — Emc — Egf) between Me and N doped GF plotted in Figj^b) suggests consistent increase in the level 
of activation of C atoms with their increasing coordination to N, which is directly corroborated by the adsorption 
energies quoted in FigJ^c). Spin-separation reduces significantly upon passivation of Co in the periodic supercell 
primarily due to restoration of full occupancy of the lone-pair orbitals of N. 

To generalize the activation of Co suggested by covalent adsorption of Me, we study adsorption of two other 
hazardous free radicals - NO and NO 2 . As evident in FigQa,b,d,e), adsorption of NO and NO 2 both happens 
through Co-N single bonds. Unlike in isolated NO, the N atom in adsorbed NO completes sub-shell filling as well as 
charge neutrality(FigQe)). In adsorbed NO 2 , both N and O atoms are charge neutral, although the N atom is over 
coordinated on account of the two N-0 double bonds suggested by the WCs, in agreement with N-0 bond lengths 
of about 1.2A. Notably, an isolated NO 2 molecule has resonating double and single N-0 bonds with the N (O) 
atom(s) positively (negatively) charged with incomplete (complete) sub-shell filling. The adsorbed configuration 
with two N-0 double bonds is stable since the energy of the electron at N in excess to what it requires to complete 
its sub-shell filling, is lowered due to its involvement in the Cq-N covalent bond. The over coordination of N is 
reflected in the long length (l.sA) of the Cq-N bond. 

Encouraged by favourable adsorption of radicals on Cq, we now probe in relation to ORR catalysis, the ad¬ 
sorption of O 2 , which is known to exist in atmosphere as diradical. EigJ^a) unambiguously suggests spontaneous 
covalent adsorption of diatomic O on Co and C’, whereas, non-nominal activation barriers appear to exist for 
adsorption of diatomic O on C atoms next to quaternary and pyridinic N atoms. Notably, pyridinic N atoms 
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at graphene edges and quaternary N atoms embedded in graphene are known to adsorb oxygen through charge 
transfer. Ground state of O 2 is a triplet where the two O atoms, each having an unpaired electron, are connected 
by an 0-0 single bond. As per WC distribution, O 2 appears (FigQf)) to be adsorb by Co through a O-Cq double 
bond.In adsorbed O 2 both the O atoms complete sub-shell filling unlike in the O 2 diradical, but both become 
charged in the process. The mono-coordinated outer O atom clearly has an extra electron (FigQf)), while the 
O atom coordinated to Co has an electron less. Thus the O 2 on the whole remain charge neutral, although their 
host Co acquires an extra electron from the rest of the GN to support its over-coordination upon adsorption of O 2 . 
The +le charge of the O atom coordinated Co justifies the longer length (1.38A) of the 0-Co double bond. The 
negatively charged outer O atom, although would attract cations, but may not be immediately able to donate its 
excess charge since such a donation would render its sub-shell filling incomplete. However, it can adsorb a radical 
and complete its sub-shell filling after donating its excess charge to a cation, or to the host GN itself which lost le 
to Co upon adsorption of O 2 . Interestingly, the 0-0 bond dissociates spontaneously upon adsorption of H, which 
can be crucial for facilitation of ORR since the outer O atom thus acts as a free active atomic O. 

Next we study the effectiveness of the experimentally observed substitution configurations to ORR through 
construction of free energy diagrams (|37l l38|) for complete reduction of an adsorbed diatomic oxygen on C’ and 
Co (Fig§. Free energy of a reduction step as a function of an applied voltage U is calculated as AG{U) = 
G(0) + neU — (|37|) . where G is the Gibbs free energy calculated as G = E ZPE — TS^nis the number 

of atomic H (proton-helectron) per O 2 : e bring the electronic charge, and ZPE the zero point energy calculated 
from the vibrational frequecies. Notably, in the free energy calculations the H atoms are considered to be at 
equilibrium with molecular H 2 . The entropy S is considered only for molecules in gas phase and has been taken 
from standard reference ([3^ . Free energy diagram for reduction of O 2 on C atoms next to isolated quaternary 
N atoms suggests (j40j) an overpotential potential of about 0.7V which is consistent with observation of ORR with 
N-doped graphene in alkaline media. However, lack of adsorption of diatomic O on such C atoms can pose a 
limitation in terms of the availability of O for ORR. C’ or Co atoms on the other hand, facilitate adsorption of 
diatomic O on those C atoms as a diradical, but tend to bind atomic O strongly enough, more than OH, to make 
the third reduction step, as described in Figj^b,c), uphill in free energy for U=0V. ORR on those C atoms thus can 
only be completed through application of a large overpotential to make all the reduction steps downhill (blue lines 
in Fig|^b, c)), which is undesirably more than the equilibrium potential of U=1.23V. Thus the thermodynamically 
limiting step for ORR on C’ or Co is the reduction of adsorbed O to OH, whereas for the C atoms next to the 
quaternary N atoms the rate limiting step is the reduction of the adsorbed OOH or OH (|4Qj) . which is also typically 
the case on metal surfaces (|38|) . Incorporating extra stabilization(j37j) of adsorbed OH due to hydrogen bonding 
provided by the water molecules in the vicinity, the overpotental requirement reduces to a value marginally less 
than the equilibrium potential. Since increase in pH from zero has same effect {G{pH) = KTln( 10)vH ^37]i . K 
being the Boltzmann constant) as that of lowering of an applied voltage (U), high overpotential implies ORR to be 
possible only in an acidic media, as has been reDorted (|4Tl [42j) . To summarize the discussion, in one hand, next to 
the quaternary N atoms, we have the C atoms which appears to support ORR in acidic as well as alkaline media, 
but may be limited by the availability of O itself due to poor adsorption of diatomic O which dissociates only upon 
adsorption followed by reduction. On the hand we have the C’ and Co atoms, which are coordinated to more than 
one N atoms, and ensure availability of O by supporting spontaneous adsorption of diatomic O, but having high 
overpotential, which stipulates an acidic media. We recall here that, experimentally, N atoms have been observed 
to substitute at next nearest sites which support existence C’ and Co, which has been also found to be energetically 
favourable. Thus it appears to be a balance between the availability of O facilitated through adequate adsorption 
of O 2 , and the ease in their complete reduction in terms of lowering of free energy, is what would finally determine 
the effectiveness of N-doped graphene as a good metal-free alternate for ORR catalysis. 

In relation to mutual covalent adsorption of N doped GNs, we now probe whether Co can be passivated by a 
similar C atom in another N doped GN. As already shown, sp^-hybridization of Co due to its passivation, would 
preserve planarity of the rest of the GF. Steric repulsion between two such planar parallel GF segments is expected 
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to make their mutual adsorption through inter-GF Cq-Co covalent bonds difficult due to proximity. However, 
additional interactions such as the dispersion interaction(D) and it — it stacking interaction are expected to help 
in inter-GF cohesion. Strength of these interactions dependent on the nature of stacking of the segments. We 
study cohesion between two isolated segments starting from an initial configuration where they are 3 Aapart. 
Steric repulsion as well as tt — tt stacking interaction both are expected to intensify with AA stacking. Interestingly, 
covalent adsorption between N-doped GF segments due to inter-GF Gq-Co covalent bond (FigJ^g)) is observed only 
with AA stacking, even without the support of dispersion interaction. Thus with AA stacking the inter-segment 
IT — IT interaction, as represented by the WF shown in Fig|^h), is stronger than the inter-segment steric repulsion. 
No inter-GF adsorption with AB stacking implies sharper drop in tt — tt interaction than that in steric repulsion 
compared to A A stacking, such that, the dispersion interaction (|3T]) alone can not overcome the potential barrier 
induced by steric repulsion. Notably, as Go can not be passivated by an ordinary graphenic C atom, N doping can 
be used to define regions of chemisorption between GF segments. 

However, with reduced steric repulsion, cohesion between N doped GNTs appears to be easier and generic. 
As evident in Figg i), unlike in case of planar segments, additional inter-GNT C-G bonds are formed in the 
neighbourhood of the inter-GNT Gq-Co bonds, implying enhanced activation of C atoms in the neighbourhood of 
Co in CNT. This is straightforward since the contractile strain caused by the Cq-N bonds of higher order causes 
reduction in surface area, which in turn increases local curvature, and thereby, the level of sp^-hybridization as 
well, leading to an enhanced level of chemical activation of C atoms in the neighbourhood of substitution by N in 
CNT than that in GF. 

To conclude, we reveal from first principles a general scenario of non-trivial chemical activation of N coordinated 
graphenic C atoms upon substitution at next nearest sites as experimentally observed. Mechanical strain due to 
increased order of C-N bonds and loss of charge neutrality of N and C atoms, in attempt to maximally retain 
TT-conjugation among C atoms, is central to the energetics of substitution by N, as well as the mechanisms leading 
to activation of N coordinated C atoms. Substitution at nearest neighbouring sites is found to be unfavourable on 
account of high mechanical stresses stemmed at the non-co-planarity preferred by the N atoms amidst strongly co¬ 
planer C atoms due to complete retention of 7r-conjugation among them. Activation is most pronounced for G atoms 
with all three nearest neighbours substituted by N possible upon tri-substitution at mutually N-Nn sites. Proximity 
to zigzag edges and curvature naturally enhances the level of activation of N coordinated G atoms. Relevant to 
oxygen reduction reactions (ORR), adsorption of diradical O 2 on activated G atoms coordinated to more than 
one N atoms, renders the mono-coordinated outer O atom reactive, although the overpotential requirement for 
complete reduction of an adsorbed diatomic O is high due to strong binding of O on such G atoms and should 
allow only an acidic media. Between two adjacent GFs or GNTs, mutual covalent adsorption exclusively among 
the N coordinated G atoms and possibly few more G atoms in their close vicinity, promises a new route to robust 
and ordered self-assembly of GNs. 

Galculations were performed partly using the central facility for high performance computing at NISER sup¬ 
ported by the Dept, of Atomic Energy of the Govt, of India (GOI), and partly in another facility supported by a 
financial grant (SR/NM/NS-1026/2011) from the Dept, of Sci. and Tech, of GOI. 
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